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ABSTRACT 

We develop the formalism for computing the magnetic field within an axisymmet- 
ric neutron star with a strong Type II superconductor core surrounded by a normal 
conductor. The formalism takes full account of the constraints imposed by hydrostatic 
equilibrium with a barotropic equation of state. A characteristic of this problem is 
that the currents and fields need to be determined simultaneously and self-consistently. 
Within the core, the strong Type II limit B <^ H allows us to compute the shapes 
of individual field lines. We specialize to purely poloidal magnetic fields that are per- 
pendicular to the equator, and develop the "most dipolar case" in which field lines are 
vertical at the outer radius of the core, which leads to a magnetic field at the stellar 
surface that is as close to a dipole as possible. We demonstrate that although field lines 
from the core may only penetrate a short distance into the normal shell, boundary con- 
I ditions at the inner radius of the normal shell control the field strength on the surface. 

Q Remarkably, we find that for a Newtonian N = 1 polytrope, the surface dipole field 

strength is .Bgurf — Hh^h/^ where Hi, is the magnetic field strength at the outer boundary 
of the Type II core and e^R is the thickness of the normal shell. For reasonable models, 
i^b K, 10^^ G and 0.1 so the surface field strength is -Bgurf — 3 x 10^^ G, comparable 
CN to the field strengths of many radiopulsars. In general, Hb and tb are both determined 

^ by the equation of state of nuclear matter and by the mass of the neutron star, but 

-Ssurf ~ 10^^ G is probably a robust result for the "most dipolar" case. We speculate 

00 

ly-j on how the wide range of neutron star surface fields might arise in situations with less 

restrictions on the internal field configuration. We show that quadrupolar distortions 
^— I are ~ — 10~^(-fffe/10^^ G)^ and arise primarily in the normal shell for B <^ Hi,. 

Subject headings: stars: neutron; magnetic fields 
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1. Introduction 

Theoretical arguments predict that the protons in the core of a neutron star form a superconductor 
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(Baym et al. 1969), which will be Type II provided that k = l/k^^p > l/v2, where S,p is the 
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coherence length and l/k^ is the London length (Easson &; Pethick 1977); numerically 
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where the gap energy in MeV is Ap(MeV) 
(because the pairing is ^5o) and Up = kpp/Sir'^ is the proton density. We expect kp^p ~ 50 
MeV at or just below nuclear density, nnuc = 0.16 F~^, which is near the boundary between the 
superconducting core and normal conducting crust inside a neutron star ( Baym et al.|197l||HaenseT 
2001) Although calculations of Ap are still being refined (Zuo et al. 



2008 



2004[ ), Eq. ^ is likely 



to be consistent with Type II superconductivity for Uh < 2nnuc- 

Two recent observations are relevant to this question. First, X-ray observations of the Cassiopeia 



A neutron star are consistent with cooling via the Cooper pair formation process (Flowers et al 



1976 Shternin et al. 2011 Page et al. 2011), and require proton superconductivity to suppress 
the contribution from URCA processes, particularly Modified URCA; if Tc^p^g ^2 — 3 throughout 



the Cassiopeia A neutron star, these processes are unimportant (Shternin et al. 2011). Second 



the discovery of a 2Mq neutron star favors relatively stiff equations of state of nuclear matter, in 
which case internal densities may not be far above nnuc (Demorest et al. 2010). We can quantify 



using relativistic polytropes, where the relationship between pressure P and baryon density rib 



is P (X n 



l+l/Af. 



to fix the scale we use the chemical potential is /in 



mb{l + enuc) at nuclear 



density, where ruh is the baryon mass. Requiring (i) 



0.065, a reasonable value (Hebeler 



et al.||2007[ ), (ii) M^ax > '2Mq and (iii) sound speed < 1 implies 0.54 < < 0.73, for which 
Mraa,x ~ 5.44(enuc)'^''^ Af© — 2.0 — 2.6M0, and central densities are n^/nnuc — 1-7 — 2.4 and 
n-b/nnuc — 2.4 — 6.7 for M = 1.4Mq and M = 2.0Mq, respectively. Actual equations of state may 
be stiffer than = 1 in some density ranges, and softer in others, but central densities may not 
be much larger than ~ 2nnuc for "typical" neutron stars. In this relatively limited density range. 



and with Tc^p^g > 2 — 3, k > l/v2 is likely to be true, so the superconductor is Type II. 

Inside a Type II superconductor, magnetic flux is concentrated within flux tubes, each of which 
carries magnetic flux $o = vr/ic/e w 2.06 x 10~^ Gcm^. The spacing between flux tubes is approxi- 
mately £b = {^o/t^B)-^/'^ = 25605^^^ F, where B = 10^2^12 G is the magnitude of the magnetic 
induction in the core; consequently 
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When ki^B ^ the magnetic field is confined to flux tubes that do not interact with one another 
to a first approximation; the local magnetic field near a fiux tube decays exponentially with a scale 
length 1/kp. Under these circumstances, the magnetic free energy is nearly proportional to B and 
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the magnetic field strength inside the superconductor is approximately ( Easson &: Pethick||1977 ) 



47r S-Kmpch? \ ^ J V50MeV/ ^ ' 

which is a function of density, since pF,p and n vary with density according to the neutron star 



equation of state; from the calculations in Zuo et al. (2008), we conclude that a fairly good approxi- 
mation is H cc with b ~ 1.6— 1.8. Since ki'^B is large, the superconducting core of a neutron star 
is in the "strong Type 11" regime. That the magnetic field strength in this regime is fixed by the 
density both simplifies and complicates the computation of magnetic structure. The strong Type II 
limit fails within a (presumably thin) transition region within which protons cluster progressively 
into nuclei until the superconducting free proton fluid disappears entirely. The outer shell of the 
neutron star is then a normal conductor. 

The challenge is to compute the magnetic field structure including both Type II inner core and 
normal outer shell, which then matches to the (vacuum) exterior. In this paper, we focus on the 
"most dipolar" external fields within stars we model as = 1 Newtonian polytropes. In order 
to find these most dipolar configurations, we require that field lines exit the core vertically. We 
then show that if Hb ^ B is the magnetic field strength at the outer edge of the Type II core, 
then the magnetic field at the base of the normal shell is very nearly —9Hb sin6, provided that the 
transition region between the core and normal shell is thin enough. This boundary condition is 
consistent with a particular dipole field solution within the normal shell. From this solution we find 
that the characteristic magnetic dipole field strength at the stellar surface under these conditions 
is approximately 

^ - ^"^^ 2.91 X lO^^G f (^(f^Y (^) . (4) 
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where the thickness of the crust is e^R and /i is the dipole moment. 

Eq. ^ is one of the central results of this paper: it relates the dipole magnetic field strength 
directly to the magnetic field in the Type II core and the thickness of the stellar crust, both of 
which may be computed given the nuclear equation of state. Although Eq. depends sensitively 
on pf^p at the boundary of the Type II core, it is noteworthy that the implied field strength in this 
model is comparable to S ~ 10^^ G, which is characteristic of many pulsars and accreting neutron 
stars. Detailed results will differ among equations of state; in particular, eb oc p^^, where pb is the 
density at the base of the normal shell, for a polytrope of index N. 

In ^we present a general overview of the problem we address here. In this section, we derive (i) 
the implications of hydrostatic balance for magnetic field configurations, (ii) the equations that we 
solve in the Type II core, (iii) the equations that hold in the normal shell, and (iv) the formalism 
for computing stellar distortions. We give specific results relevant to the "most dipolar" case and 
derive, in particular, Eq. (Q, but also highlight causes of perturbations around this simple solution. 
In ^we illustrate some of the salient features of the problem via a toy model that is totally analytic 
and surprisingly close to being realistic. Then in ^we present particular results for the most dipolar 
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case of greatest interest. In this section, we consider not only the simplest version of the "most 
dipolar" case but also perturbations around that solution that lead to non-dipolar corrections to 
the fields. We conclude in ^by reviewing these solutions and suggesting extensions to other cases 
that may allow surface field strengths that are either appreciably lower or higher than Eq. Q. 



Overview 



In this paper, we compute the magnetic field for a star with a superconducting core that matches 
onto a normal crust surrounded by vacuum. We model the star as a Newtonian = 1 polytrope 
with is a poloidal magnetic field that only distorts the star slightly. Previous calculations of this type 
were done for poloidal fields of a uniform density star ( Roberts|1981 ) and for toroidal magnetic fields 
in a Newtonian A'^ = 1 polytrope (Akgiin & Wasserman 2008, Akgiin 2007). These two different 



previous calculations were simpler than those presented here. A significant complication in finding 
the magnetic field structure is that H is density-dependent according to Eq. accommodating 
this complication will require a generalization of the method of solution for a uniform density star 
(Roberts 1981). Calculations for a toroidal field could account for the density dependence of H 



because it was unnecessary to determine the field line shapes, which were specified a priori (Akgiin 



& Wasserman 2008) 



Our calculations are perturbative in that we assume that the distortions induced by magnetic 
stresses (which we compute) are small. However, we take full account of the requirements of 
hydrostatic balance, which constrain the field shapes, just as for normal conductors ( Wentzel||1961 



Prendergast 1956); our calculations are analogous to those done previously for small magnetic 



distortions in a normal conductor ( Monaghan|| 1965 ). We find the "true equilibrium" configurations 
in which the magnetic distortions are assumed to obey the same equation of state as the unperturbed 



star. The actual magnetic fields in neutron stars may not be true equilibria in this sense (Mastrano 



et al. 2011) even though, on sufficiently long time scales, the fields should relax to these states 



naturally. Even in non-equilibrium field configurations the field at the base of the normal shell will 
have to match properly to the field in the Type II core. We therefore expect a version ofEq. (|4j) 
to remain true for substantially dipolar fields even for configurations that are not true equilibria. 

We focus on the "most dipolar" field configurations that arise when field lines hit the outer edge 
of the Type II core vertically. As we shall see, this is a special configuration, and infinitely many 
others are possible. 

Throughout this paper, we assume that the entire core of the star is superconducting. We argued 
above that this may be reasonable if the density range in the core is not too large. We expect 
that the external dipole field that emerges from our "most dipolar" solutions would not be affected 
significantly if the very inner core of the neutron star is not superconducting. 
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2.1. Requirements of Hydrostatic Balance 



The magnetic free energy is /(p, B), where p is mass density and B is magnetic induction; then 
the magnetic field is H = 47r9//5B, and the magnetic force density is 



^mag 



dl 
dp 



+ 



(VxH)xB 



-pV 



dl 

dp 



+ 



JxB 



(5) 



using Ampere's law, VxH = AirJ/c (Mastrano et al. 2011). The equation of hydrostatic balance 
including pressure, gravity and magnetic forces is 

V p fmag 

= -—-V^+ (6) 

P P 

where ^ is the gravitational potential, P is the pressure and p is the mass density. For a barotropic 
fluid P = P{p) and dh{p) = dP{p)/p, in which case 







V /i + ^' + 



dl 
dp 



+ 



JxB 

pc 



(7) 



Eq. ([T]) is only consistent mathematically if JxB/pc 
axisymmetry, 



B = Vx 



r sin( 



and H 



Hp + Ht4>] the current density is 

VxH = VxHp + -^(^-^^^^^)^'^ 



— V<^, where $ is some potential. In 

+ BT{r,9)^ (8) 

_ 47rJ(^ V(-ffT''sin6')x<^ 



+ ' ^ . n' ' ■ (9) 
r sm tier sm 

Requiring that </)'(JxB) = implies that HttsvuO = T-LxiA); the Lorentz acceleration JxB/pc is 
a total gradient if 

cBt d-HriA) 



J 



A-K dA 



cpr sm9J'{A) 



(10) 



Eq. (10) is well-known for normal magnetic equilibria (Prendergast 1956 Woltjer 1960 Wentzel 



1961t |Monaghan||1965| ); the Bt = version of Eq. alOh has also been derived previously for Type II 



superconductors with poloidal fields ( Roberts|1981 Akgiin|2007 ) . Similar results have been derived 
independently by Lander (2012). 



We are interested in poloidal fields only, so Bt = 0. (See also Lander 2012, for specific models 
including poloidal and toroidal fields.) We also assume that /(/O, B) = f{p,B) i.e. there are no 
preferred directions in space so the free energy only depends on i? = |B|. In this situation, the 
magnetic field and magnetic induction are parallel: H = AirBdf /dB where B = B/B. 

Our goal is to solve 

VxH = Vx{HB) = VHxB + HVxB = Anpr sin 9 J (A) 4> (11) 

throughout the star. We shall strive for solutions without surface currents; in particular, we are 
interested in solutions in which field lines cross the equator smoothly and vertically. With this goal 
in mind, we will have to determine both the field configuration and its source J{A) self consistently, 
within both the Type II core and the normal shell. 
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2.2. Field Configuration in the Type II Core 



In the strong Type II limit, H = H{p) in the core, and invoking the assumption that magnetic 
distortions are small, we can substitute the unperturbed density of the background star, /o(r), so 
H = H{r). 



With this simplification, Eq. (11) is a complicated partial differential equation even if we specify 
J{A) somehow, and we shall not attempt to solve it directly. Instead, we extend the procedure 
first employed by Roberts (1981) for uniform density stars to stars with H(r) ^ constant. Since 
Bt = 0, Eq. A8n implies that B-^A = i.e. A is constant along poloidal field lines. This allows us 



to consider Eq. (11) along individual field lines. 

Let us concentrate on a particular field line. Introduce the parametric independent variable s, 
the arc length along this field line; position along the field line is r(s) and the tangent to the 
field line is B(s) = dr{s)/ds. There are two directions normal to the field line: cf) and a second 
direction h{s) which we define to be n{s) = (f)XB{s). The field line curvature is K{s) given by 
K{s)h{s) = (f'r{s)/ds^ = B{s)-'VB{s). Eq. (11) is then equivalent to 



AttJ 



A-npr suiej{A) 



-VH-h + KH . 



(12) 



For uniform H the first term in Eq. (12) is absent and we recover the equation found by Roberts 



(1981). In Eq. (12) J{A) is simply a constant associated with the field line. 

In axisymmetry, we may choose to specify position along a field line by r(s) and 9{s); in that case 



Bis) 



. dr(s) 
ds 



+ Oris 



dejs) 

ds 



r cos A(s) + sin A(s) 



and therefore 



n{s) = 4»<B{s) = -fsinA(s) + 0cosA(s) 



The field line curvature is 



K{s) 



dA sin A 

ds r 



using Eq. (14) and Eq. (15) in Eq. (12) we find 



dA 

ds 



Anpr sin 9 J' (A) sin A / r dH{r 



H{r) 



H dr 



+ 1 



(13) 

(14) 
(15) 

(16) 



Eq. (16) and Eq. (13) determine the field line shape jointly, given p{r), H{p) and a value of J {A). 

For our solutions we shall assume an = 1 polytrope for which p[r) = p(0) sinx/x, with x = nr / R] 
we shall also assume H oc p^ . With these choices, and the definition ds = Rda/Ti, Eq. (16) becomes 



dA 
da 



where 



J'(^)x''sin6' sin A 
(sinx)^^^ X 



[6(xcotx - 1) + 1] 



m^p{0)J{A) 
ttH{0) 



(17) 



(18) 
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the field line shape is found by solving Eq. (17) along with 

dx d9 

— = cos A X— = smA . (19) 
da da 

We are interested in solutions that start at some equatorial footpoint where x = x(0) and 9{0) = 
7r/2. At this footpoint, symmetry dictates that the field line be perpendicular to the equator; 
choosing the direction to be vertically upward implies A(0) = —tt/2. 

For the "most dipolar" case we also want field lines to hit the outer boundary of the Type II 
region at x = x^ with A = —9. Very generally, it is impossible to find solutions that are exactly 
vertical throughout the core for a given H(p). (The solutions are nevertheless close to the toy 
model developed in ^that demands that field lines are exactly vertical.) The condition that a field 
line starts out at x(0) and 9 = ^/2 pointing vertically and also hits vertically at Xf, determines the 
value of J {A) for that field line. Note that in this context, is just a label for the field line; we 
could just as well label the field line by x{0). In fact, this procedure does not determine A{x,9) 
within the Type II core. The end result is J'{x{0)) and, since field lines hit Xb at 9{x{0)), Jb{9) is 
determined parametrically on the outer boundary of the core. 

We can also consider other types of solution for field lines, for instance solutions in which field lines 
hit Xh with some other specified set of orientations and, therefore, carry different values of J{x{{))) 
along with them. We must be careful to choose orientations for which field lines do not cross within 
the Type II core. But even that restriction allows many different possibilities other than our "most 
dipolar" case. 

The fact that we do not determine A{x, 9) in the Type II core means that even with a given set 
of solutions for individual field lines, we can consider various "field line densities" within the core. 
This freedom can accommodate theories of field line evolution in which secular motion of flux tubes 
leads to different concentrations of flux within the core. 



2.3. Solution in the Normal Shell 



Within the normal shell, it is most straightforward to solve Eq. (11) directly: using H = B and 
substituting for the poloidal field we get 



F(0)i?2 



[J'{A)x sin x sin^ 9] 



d'^A 1 fd'^A 
+ - 



di 



89^ 



cot 9 



dA 
~d9 



(20) 



where we have used the same parametrization as in Eq. (18) and specialized to the = 1 polytrope 



density profile. Eq. (20) would be straightforward to solve given J {A). However, we have to 



determine J {A) to be consistent with matching conditions at both the inner boundary of the 
normal shell at Xh and at the stellar surface, where the field must match smoothly to vacuum. 

Matching to a vacuum exterior is accomplished most easily by introducing the expansion 



A{x,9) =s\n9^Aj{x)P2j+i^{cos9) 

3=0 



(21) 
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where P2j+i^{p) is an associated Legendre polynomial. Using Eq. (21) in Eq. (20) implies 



d^Ajix) {2j + 2){2j + l)Aj{2 



+1 



d^l^/l^P2J+l\^Ji)J{A) 



where 



2(2i + 2)(2j + l) 
4j + 3 



f22) 



(23) 



The source term in Eq. (22) generally involves the entire set of Ai, not just i = j. But once J{A) 
has been determined, the source term is known, and the boundary condition at the surface x = ir 
is 

(24) 



dx 



for matching to a vacuum field in the exterior. 



dA,{x)\ (2j + l)^,(7r) 



TT 



Matching to the core requires that the normal component of B and tangential component of H 
be continuous across the boundary. In the simplest case, the transition zone thickness may be 
neglected. If Ab{6) is the value of A for a field line hitting Xb at 9 and B{6) is the magnitude of the 
core induction field there, then the matching conditions are {/j, = cos 6) 



fdA 

dA\ 
9x / 



R^xlB (9) cos Abie) 



sin sin Ab{9) 

^2 



(25) 



In the strong Type II regime, B{9) <^ Hb and unless Ab{9) < B/Hb the field is predominantly 
tangential at the base of the normal shell. To a first approximation, {dA/dfi)b ~ 0, which is 
equivalent to Aj{xb) ~ 0, in the strong Type II regime. 



Eq. (24) and Eq. (25) impose three conditions on each Aj{x), which is one too many. This overde- 
termination implies constraints on the form of J^{A) within the shell. 



In the "most dipolar" case, Ab{9) = —sm9, and Eq. ( |25[ ) implies that the field at the inner edge 
of the normal shell is approximately —6Hbsm9 in the limit i? — t- 0. In this limit, field lines from 
the interior do not penetrate into the normal shell but still infiuence conditions there via Eq. (25). 
With B = 0, the field inside the normal shell is precisely dipolar, A{x, 9) = Aq{x) svc? 9, where 

HbR^f{x;xb) 

(Sx^ — 6) sin X — {3x1 — 6) sin Xb 



Ao{x) 
f{x;xb) 



3^2 

4 



X? + 6x cos X — 6xb cos Xb 



x{l + COSXfc) 



(26) 



Within the shell, consistency with Eq. (24) and Eq. (25) require that 

(sinxb/xfc)* 



J (A) = J^shcU 



Hb 



i7(0)(l + cos Xft) 1 + cosXb 



(27) 
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which the simplest case showing how the boundary conditions determine J {A). The steUar dipole 
moment is 



if vr 



Xb 



fi = RAo{7r) -- 
neb < 1, /(vr;xb)/37r2 ^ Sb/Sn 



37r2 

which leads to Eq. Q for thin shells. 



(28) 



Eq. (26) holds when we can neglect the thickness of the transition layer and the radial field that 



pokes into the normal shell through it. If the transition layer is thin enough, then its effect is to 
introduce a "surface current" into the jump condition from the Type II core to the normal shell: 
the tangential field at the base of the normal shell is now 



Be{xb, 



vr 



Rh 



"Xb smf \dx J 
where i is the thickness of the transition zone. 



Hb sin Abi9) + 



Trejb{e)H{0)smxb sin (9 
R 



(29) 



1 contains a brief discussion of the rather complex 



physics of this region. There, it is shown that Eq. (29) holds as long as field lines do not rotate very 



much as they pass through the transition zone. Roughly speaking, this will be true provided that 
Hbi/BR <C 1; therefore the second term in Eq. (29) is ^ BH{0)/H^ times the first. Nevetherless, 



Eq. (29) perturbs the solution away from Eq. (26) even in the "most dipolar" case. In particular 



the jump term engenders non-dipolar fields in the shell and outside, and requires a more detailed 



calculation of J^{A) inside the shell. As long as the corrections to Eq. (26) are small they may be 



computed perturbatively. A method of calculation is outlined in Appendix |Xj 

Relaxing the i? — J- limit also perturbs the solution even if the transition zone thickness is negligible. 
If we assume that 



A{xb, e) = \Brl sin2 9 = 

then the unperturbed field inside the normal shell remains dipolar, and is altered to Aq{x, 
Aq{x\ B) siv? 9, where 



Ao{x;B) 



BR^xl ^ {Hb + B/2)R^f{x;Xb) 



27r2 



X 



37r2 



and the current density parameter is 

c7shcll 



Hb{l + B/2Hb) 
F(0)(l + cosxfe) 



Eq. (31) implies a magnetic moment 



Brl , {Hb + B/2)R^f{TT-Xb) 



+ 



3^2 



(30) 



(31) 



(32) 



(33) 



which reduces to Eq. (28) for = 0; the first term in Eq. (33) is simply the dipole moment 



associated with B and the second arises from the currents in the normal shell. 



Eq. (|31|) implies that magnetic field lines from the core enter the shell. The bounding field line 

^b^ 



entering from the core has A 



^Br^, and therefore defines a region 9 < 9u{x) where 



sm 



hiix) 



/3 



f3xb/x + f{x;xb) 



(34) 



9 



in Eq. (pll) 



a _ ^^^Brl ~ _ f{x;xb) , . 

The region occupied by the field hnes is small if /3 ~ B/Hbeb ^ 1; note that this is more restrictive 
than the condition for validity of the strong Type II regime, B <^ Hb- Eq. ( |34[ ) indicates that the 
bounding field line rises from 6 = t: /2 close to Xb, remaining at ~ 1 as long as f{x;xb) ~ /?, which 
is the case for x/xf, — 1 < /?. Ultimately, the bounding field lines emerge at x = vr confined to a 
polar cap with 6 < The volume of the region occupied by impinging field lines is only ~ /3 
times the volume of the shell. In the limit /? — )• field lines from the core do not penetrate into the 
normal shell. 

However, within that volume, AJ'{A) = J{A) — J^heW 7^ because field lines entering from the 



core carry their j{A) along with them. The zeroth order solution represented by Eq. (31) and 



Eq. (32) requires perturbative corrections driven by the difference /^.^{A) within this region; the 
corrections are 0{I3/S.J{A)/ Jshcw)- Appendix [a| outlines how these corrections may be taken into 



account. 



2.4. Magnetic Distortions Due to Poloidal Fields 

For poloidal fields, the Lorentz acceleration is — V<I> = i7(A)VA, so 

7rF(0) '-^ 



HA) 



dA'j{A') 



and the Bernoulli equation is 



Of 



C = h + '^ + ^ + 
op 



(36) 



(37) 



where C is the Bernoulli constant. For perturbations of a = 1 polytrope, if we expand in 
Legendre polynomials 

^ix,e) = Y,^iix)Piicos0) (38) 



we get 



LA. 

x^ dx 



dx 



+ 



^!i{x) = 5C5i,a 



dl 
dp 



+ $ =SC5ifi + Miix) 



(39) 



where the magnetic potentials are expanded in the same fashion as Eq. (38) and 5C is the pertur- 
bation of the Bernoulli constant away from its background value. The general solution of Eq. (39) 
that is regular at a; = is 

nx 



^i{x) = 6C5ifi+ji{x) 



Ki 



dx'{x'fyi{x')Mi{x') 



+ yi{x) / dx'{x'fji{x')Mi{x') (40) 







where ji{x) and yi{x) are spherical Bessel functions. The constants Kq and 5C are determined by 
the conditions = = [d^!{x)/ dx]-,^^ which imply 

1 



5C = Kq 



dxx jo{x)M(j{x) 



(41) 
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and the constants Ki are determined by the conditions [d'^i{x)/dx].jr + {I + l)'^i{ir)/'!r = 0, which 
imply (//(x) = dfi{x)/dx) 



= Ki 



+ 



(^ + i)yKvr) 



TT 



dxx^ji{x)Mi{x) . 



For I > 0, 



= 2-/r ■ / ^ / dxx^jiix)Mi 



(x). 



(42) 



(43) 



We may define the stellar multipolc moment by ^/(vr) = —GMAii/R. With this definition, the 
contribution to the moment of inertia tensor resulting from the quadrupolar distortion M2 is 
5I2 = -^Mi?2>f2(-xx - yy + 2zz). 

In the "most dipolar" case with negligible transition layer thickness and core magnetic induction, 
the distortion arises entirely from magnetic stresses inside the crust and is quadrupolar. In this 
case, (6(^;) = 1 for 2; > 0, zero otherwise) 



M2(x) = - 



187rp(0)(l + cosxfe) 97r2M(l + cosa;b) 



(44) 



and (tt ^2 (tt) + Stt ^'2 (tt) = vr) 



2/72/^4 dxx''j2{x)f{x;xb) 



r dxx^'j 

4 1 



+ cos X}, 



2H^R^I{xb) 



(45) 



For thin shells, T{xb) ~ 97r/2 = 14.14, independent of the shell thickness; for = 0.1, it is 
X(0.97r) = 11.44; consequently 



7W2 = -1.39 X 10-^ Hli^RtoM^i 



I{xbl 
10 



(46) 



where Hb = 10^^ Hb^u G, R = 10i?io km and M = IAMi^^Mq. Thus, the scale of the quadrupolar 
distortion is determined by the value of Hb- 



3. Toy Model: Vertical Field in the Core 

We shall be interested in solutions that lead to external fields as close to dipolar as possible; such 
solutions have fields that are vertical at the outer boundary of the Type II core. In general, these 
will not allow field lines that are exactly vertical throughout the core. 

Nevertheless, to develop a toy model, suppose H = H{r)z; this model will illustrate many features 
of the problem. For vertical fields, Ju{A) = constant = Ju and from Ampere's law 

TT 
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where x = irr/R and H{R) = 0. The dimensionless combination 

4 i?V(0)J ii _ , 

moreover 



Jii ^ " = \ ; (48) 



(ilnff xsmx 

6 = — = -. (49) 

amp (1 + cosa;)(l — X cot x) 

varies monotonicahy between 6 = 3/2 and 6 = 2 from the center of the star to the surface. Although 

this model is not truly physical because H[p) = H{p/p{0)) depends on the central density, hence 

mass, of the star and so involves a density scale that is not the same for all stars, 3/2 < 6 < 2 is an 

acceptable range for nuclear equations of state. Moreover, when we consider models with H (x 

with realistic values of 6, we shall see that even though J'{A) is no longer constant, its values are 

close to 0.5 for cases where field lines hit the outer boundary of the core vertically; because field 

lines are vertical at the equator, they remain close to vertical throughout the core in these cases, 

as may be seen from Fig. [T] 

The interior field must match onto the field in the normal shell. If the transition from supercon- 
ductor to normal conductor occurs in a transition region of negligible thickness, then the field is 
dipolar in the shell to a good approximation, and is given by Eq. (31) and Eq. (32). Eq. (32) holds 
generally for the dipole field within the normal shell; for the toy model Hb = ^H{0){1 + cosx;,) so 

j.«> = j(l + 5|^)=j„(l + ^) . (50) 

The core and shell current density parameters are nearly the same in this toy model for B/2Hiy <^ 1. 
This will also be true for more realistic models, but the difference will be ~ 0.1, not ~ B/Hf,. For the 
toy model, the perturbations caused by AJ'{A) are very small, since (3{Jii — jfsheii) = \PBJi\/Hi,, 
and we shall ignore them; for realistic models, where Jn — Js'hcW ~ 0.1, the perturbations are larger, 
~ 0.1/3. 

In general, there will be a contribution to the stellar distortion oc B. One reason is that H^, — ^ 
i?fe + i?/2 in the dipole solution within the shell; thus we substitute —?- {Hf, + 3/2)"^ ~ H^ + HfyB 
in Eq. (45). There are also contributions from the terms in the solution that are oc B. For the toy 
model with = ^ — i^sheii the extra / = 2 contribution to the driving term is 

24/t^/9(0j \ Xb/X it Xb < X < IT 

which implies an additional quadrupole moment 

H{0)BrlR^AI{xb) 1-03 x Hl,,Ri,AIixb) B 

^^' = 6^GM^ = Ml. ^ ^ ^ 

where AX{xb) = AXii{xb) + AXsheii(2;fe) with core and shell contributions 



2{rb/Rf 

AXii[xb) - 



1 + cos Xb 



15 \ / 15 , 

^ — 6 sm Xb+ \ Xb cos Xb 

xt J \ Xb 



AX^Yic\\{xb) = ^("^^Z"^-* [-Xb{l + cos Xb) + 3 sin Xb] (53) 

1 + cos Xb 
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respectively. The crust contribution AXsheii(a^fe) is the same in reahstic models. For thin shells, 
AXsheii{xb) — 12/5;, and AXsheii(0.97r) = 26.1. The core contribution AXn(xfe) will be close to 
the value in the toy model, with differences arising because field lines are not strictly vertical, 
realistically. For thin shells, AIii{xh) — {60 / tt — Air) / 6^ ~ in the toy model and AXn(0.97r) = 

35.8. Comparing Eq. (52) with Eq. (46) we see that |AA^2| ^ {(3/eb)\M-2\ characteristically. 



4. Numerical Solutions 
4.1. Field Line Structure 



In the strong Type II limit where H = Hci{p) does not depend on B, we can solve Eq. (17) and 



Eq. ( 19 ) for each magnetic field line individually. Symmetry dictates that field lines be vertical at 
their equatorial footpoints. For the "most dipolar" solutions that match to a pure dipole field in 
the normal shell to zeroth order in i and we also require that each field line be vertical at Xf,. 
We can label each field line by its footpoint radius, xq. Once we specify H{p) = Hci{p) we can find 
the value of j{xo) for which these conditions are satisfied. Thus, we determine both the source 
and the field line shapes that are consistent with our boundary conditions. 



In order to choose the function H{p) we fit the values of the proton fraction given in Zuo et al. 



(2008) to power laws of the baryon density. Writing H{p) oc p'^, we found that 6 1.6 — 1.8. We 
used the A'^ = 1 polytrope for the mass density profile of the star. 

Fig. [TJshows the results for b = 1.6, 1.7 and 1.8. Note that there are two competing features of the 
solutions: Vi^ tends to curve the field lines inward, and J^{xo) opposes this trend. The field lines 
are most nearly vertical throughout the core for b = 1.6, the smallest value displayed. Even though 
the field line shapes are not precisely vertical throughout the core, the deviations are relatively 
small. We have also found the field line solutions for 6 = 1, 1.5, and 2; they are substantially 
similar to the configurations shown in Fig. [T} 

The lower right panel shows the current density parameter J/j, as a function of sin^ 0. This quantity 
is found by determining j{xQ) for each field line; using the solution for 0{a) along the line, we find 
where it intersects the surface x^. Note that our solutions do not determine the value of A for a 
given field line. 

For B = = P and £ = boundary conditions at the core-shell interface dictate that B = 
—OHijSmO at the base of the shell. This means that in this approximation field lines do not 
penetrate from the core to the shell. The solution inside the shell that matches smoothly to a 



vacuum field is given by Eq. (26). Since there are two boundary conditions at Xb, the surface 



boundary condition determines Jsheii (Eq. (27)) so the entire solution outside the core not only is 
determined fully but also depends on Hb- 

Fig. [2] illustrates the field line configuration for 6 = 1.6 but (3 = 0.1, to lowest order in f3. 



(Corrections will be considered in ^4.3 ) In this case, field lines from the core penetrate into the 



normal shell. The inner thick line shows the bounding field line. It is apparent that field lines hug 
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Fig. 1. — Field line configurations. Dashed lines at x = O.Qtt and x = it indicate the outer boundary 
of the Type II core and the stellar surface, respectively. The lower right panel shows Ji, = J{xq) 
as a function of sin^ ^ on a; = x^. 
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close to X = Xh before fanning out into a polar cap whose angular extent is ~ v73. Also shown in 
this figure is the small region where field lines in the shell are closed: the first open field line has 
its inner footpoint at a; ~ O.OTvr. (This region is also apparent in the /? = configurations in Fig. 
[Tj) In a fiuid shell, this region would be magnetohydrodynamically unstable ( Wright || 19 73 Markey 
I&: Tayler 1973) but shear stresses in the normal shell restore stability since the shear modulus 



lO^^iJfagcm ^s ^ (e.g. Haensel 



2001 



Horowitz &; Kadau 



2009 



Baiko 2011)). 




Fig. 2. — Field line configuration for /3 = 0.1. 



4.2. Magnetic Distortion: Mass Quadrupole Moment 

In general, there are four contributions to the quadrupolar distortion to linear order in B/Hh'. 



M2 = -l^-^HhAM^l 



Ci + {C2 + C3 + 



B 



(54) 



The coefficients Ci and C2 arise from the shell, and are independent of the core structure; for 
rft = 0.9i2, (Ci,C2) = (+1.58, +28.3). The coefficient C3 arises from ^{A) and the coefficient C4 
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arises from Bdf /dp; these depend on the structure of the core field: for = 0.9R 



(C3, C4 



(+30.4,-0.434) if 6 =1.6 

(+35.4,-2.20) if 6 =1.7 

(+40.7,-4.64) if 6 =1.8 

(+36.7,0) Toy 



(55) 



where in all cases C3 > IC4I because field lines are nearly vertical (exactly vertical in the toy model). 
Roughly speaking, the coefficient of B/Hb in Eq. ([54| is ~ 60, but B/Hb « 2/(7r; Xb)R^(3/3Tr^rl = 
/3/14.1 for rb = 0.9R, which mitigates the dependence of A^2 on B for small f3. 



4.3. Perturbative Corrections: Nonzero i and /3 



Table [T] shows the results of computing perturbations resulting from nonzero transition zone thick- 



ness using the method developed in Appendix A.l Although the table only contains results for 
multipoles Qj with j < 5, the calculations were performed using J'bj with j < 14, and therefore 
determined Qj for j < 14. The procedure can be iterated but we did not do so. To obtain these 
results, we solved Eq. (A12) with AAj{xb) = 0. Appendix [B| shows that in fact AAj{xb) 7^ 0. 



However, according to Eq. (B7), AAj{xb) is primarily dipolar to 0{£). This dipolar component as 



well as the effect of J^bfl may be absorbed into a "renormalized" background solution. The second 
line in Table [T] tabulates the "extra" change in Qq that results when we renormalize in this fashion. 



3 









-4.918e-01 


1.701e+01 








5.239e-01 


1 


6.067e-03 


-4.999e-01 


2 


-5.198e-03 


2.723e-01 


3 


1.241e-03 


-6.844e-02 


4 


-3.843e-04 


1.973e-02 


5 


8.526e-05 


-3.657e-03 



Table 1: Values of Jb,j (col. 2), Qj{i)/R^Hp (col. 3) for j < 5 (col. 1) and b = 1.6 and n = 0.9i?. 
Second line renormalizes Hb as described in the text. 



Table [T] shows that the multipole perturbations that result from nonzero £ are small but not entirely 
negligible. For example, the shift in the dipole moment is Apj [i = —ILQli/R; for l/R = 0.01 this 
is a downward change of about 17%. Since we require AAj[xb) to be relatively small, Eq. (B7) 
indicates that i/R <^ B/Hb ~ /^Cfo, so consistency would require somewhat smaller — A/i//i: shifts 
are restricted to < a few percent. The other multipole moments are smaller by at least a factor of 
30. 
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Implementing the procedure in Appendix |A.2.1| proved more challenging computationally, partic- 
ularly for larger values of /3, where the various matrices that had to be inverted turned out to be 
nearly singular. These computational problems did not arise for /3 < 0.05, and for 10~^ < /3 < 0.05 
we found A/x/|U ~ —0.115/3 and a RMS change AAj.^s{R) / Ao{R) — 0.150/3 assuming b = 1.6 and 
€b = 0.1. Thus, although the perturbations are predominantly dipolar, nondipolar contributions 
are only a few times smaller. These results were found using various maximum < jin < 10 in 
the expansion for i7b(^), and various maximum 10 < jmax ^ 14 in the multipole expansion within 
the normal shell. Empirically, we found that about 90% of the perturbation arose when we took 
jin = 0, which suggests that the main effect arises from the peculiar shape of the region within the 
shell occupied by field lines entering from the Type II core; see Fig. [2} 



5. Discussion and Conclusions 

In this paper, we developed the formalism for computing the configuration of poloidal magnetic 
fields inside a neutron star with a strongly Type II superconducting core. Qualitatively, the re- 
striction to hydrostatic balance without "magnetic buoyancy" imposes constraints on the electric 



current density given by Eq. ( 10 ). Because of these constraints, we must determine both the current 
densities and field line structure at the same time. Within the Type II core, the task is simplified 
when B <^ H, for in this limit H is a function of density, and therefore a function of radius, to 
lowest order in the magnetic distortion. 

The main results of this paper are: 

1. The Type II core provides an important boundary condition for fields in the normal shell 



surrounding it; see Eq. (25) 



2. If the core is strongly Type II, continuity of the tangential and normal components of magnetic 
field imply that the field points very nearly along the inner boundary of the normal shell, 
unless the field just inside the boundary is tuned finely. Thus, generically field lines from the 
core do not penetrate deep into the normal shell of the star. However, the field strength at 
the base of the normal shell is ~ Hi, in this case. 

3. Substantial currents in the normal shell are required for the field to adjust in strength and 
direction to conform with vacuum boundary conditions at the radius of the star. 

4. In the "most dipolar case" in which the field is vertical at the outer edge of the Type II core. 



the emerging field is dipolar, with a surface field strength given by Eq. (26). 



5. If the thickness of the normal shell of the star is €bR with ef, <^ 1, the surface field strength for 
the "most dipolar" case is ~ H^eb/'^- For typical neutron star parameters Eq. ([4]) implies a 
surface field strength ~ 3 x 10^^ G, which is characteristic of the large population of "normal 
neutron stars." 
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6. Magnetic fields induce a quadrupolar distortion of the star — — 10 ^. The scale of the 
distortion is set primarily by Hf, for small B/Hi,; see Eq. (54). 



7. Perturbations arise from the finite thickness of the transition layer in which the superconduc- 
tor disappears and from field lines that penetrate from the Type II core to the normal shell. 
These perturbations are small if i/etR ^ ^ 1, which delineates the domain of validity of 
our calculations. 



Although our calculations assume a star in hydrostatic balance, some of the qualitative features 
found here would hold for stars that are out of equilibrium as well. For example, in configurations 
without surface currents, we would still expect field lines to emerge from the core pointing largely 
along the inner boundary of the normal shell, with a field strength there ~ Hh. Consequently, we 
would still expect the emerging field to have a field strength similar to what is found in equilibrium. 

Our calculations have employed a very simple equation of state: the = 1 polytrope. Within the 
context of this equation of state, if we stipulate that pb is fixed by microphysics of the phase change 
from core fiuid with free protons to a crystalline crust with protons confined to nuclei, then Hh will 
be independent of the stellar mass. However, since R is fixed for a (Newtonian) N = 1 polytrope 
the central density is tx M and the fractional thickness of the normal shell is oc 1/M. Thus, we 
would expect some variation of the surface magnetic field with stellar mass, but this would only 
lead to a rather small range of values differing by less than a factor of two. 

More realistically, with Hh fixed the surface field strength will still be cx Hbeb, but the value of 
eb and the proportionality constant would depend on details of the equation of state. We expect 
this to remain so even if the Type II region does not encompass the entire stellar core, as we have 
assumed here: Type II superconductor with an inner boundary would have different currents but 
the lower boundary condition on the normal shell is sensitive primarily to Hb- Nevertheless, very 
generally we expect that the "most dipolar" case will result in field strengths ~ 10^^ G with little 
variation. This is contrary to observations, which show that even the dipolar field strengths, as 



indicated by the P — P diagram for radiopulsars (e.g. Lorimer & Kramer 2012), vary over several 
orders of magnitude. 

The obvious solution to this problem is that the magnetic field configurations for many pulsars are 
very similar to the most dipolar configuration, but that there may be substantial deviations. We 
envision several reasons for such variation: 



1. The strong Type II limit fails because i? ~ if in the core. We expect that this is true for 
magnet ars, for example. 

2. Failure of the strong Type II limit is an extreme example of a case that is outside the scope 
of the calculations in this paper. There are configurations with B <^ Hb that violate our 
requirements: (a) If e^^ ^ /3 > 1 field lines impinging from the core affect the surface field 
non-perturbatively. (b) If f3 is small enough that i/R > B/Hb ~ f3eb then field lines curve 
significantly within the transition region between the Type II core and normal shell, which 
would affect their coupling. 
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3. Toroidal magnetic fields, which were not considered here, can alter the poloidal field configu- 
ration even in hydrostatic equilibrium. We will treat this possibility elsewhere. Toroidal fields 
also mitigate instability. We have seen that there is a region in the normal shell that would 
be unstable if the normal shell were a pure fluid (see Fig. [2]) although in a crystalline crust 
shear stresses should suffice to ensure stability there. However, the Muzikar-Pethick-Roberts 
instability in the Type II core ( Muzikar &: Pethi"ck||1981 Roberts||1981 Akgiin &; Wasserman 



2008 ) might require toroidal fields to ensure stability. 



4. Even if none of the last three caveats holds, so the core is in the strong Type II regime and 
toroidal fields do not affect the poloidal fields significantly, the requirement that field lines hit 
the outer boundary of the Type II core that is the basis for the most dipolar configuration is 
very special. Relaxing it could lead to greater diversity in surface field strengths. 



We expect that of these options the last is most important for neutron stars with surface magnetic 
field strengths substantially below 10^^ G. The extreme case is when field lines hit the outer 
boundary of the Type II core radially. In this extreme, the field remains radial as it enters the 
normal shell, and is weak in strength, ~ S ^ Hi,. Entering field lines carry current density 
~ cpbH(0)/ p(0)R with them, and therefore curve downward toward the stellar equator within a 
short distance of the boundary: crudely the field should penetrate a distance ~ Bp{0)R/ pi,H{0), 
which is ~ I3{ph/ p{0))^~^ <^ 1 times the thickness of the shell. We might then expect that the 
entering field lines are confined to a small region of the normal shell, and the field strength at the 
surface could be < -Bet, which may be very small. We will report elsewhere on calculations that 
employ different boundary conditions at the outer radius of the Type II core. 

Even though it relies on special conditions, the "most dipolar" solution reveals that there is a 
characteristic surface magnetic fields strength that is determined by the nuclear equation of state. 
Our calculations show that this field strength is ~ reassuringly close to the values deduced 

for many neutron stars. By relaxing the restrictive boundary conditions underlying the most dipolar 
solution, it seems plausible that a variety of field strengths may be attained in equilibrium. Field 
strengths > 10^^ G or < lO^^ q 

may be attained under other restrictive conditions. 



We thank S. Lander for helpful correspondence. This research was supported in part by NSF grant 
AST-0606710, by a fellowship from the NASA/New York Space Grant Consortium and by the 
College of Arts and Sciences, Cornell University. 



A. Perturbation Solutions 

We begin with a zeroth order solution, Ao{x,6), that corresponds to a given i7shcii(^); we use Eq. 
(26). Perturbations distort the solution to 

oo 

A{x,e) = Ao{x,e) + sm9^AAj{x)P2j+iHcose) . (Al) 

j=o 
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The perturbed solution corresponds to jiA) = jshdi + (A) , where Aj'{A) must be determined 
along with the perturbed field solution. To do this, we begin by choosing a functional form for 
Aj7(j4) which will necessarily involve unknown parameters to be determined. Since A Aq hy 
assumption, we substitute AJ'{A) ^ AJ'{Aq) to find the perturbed field; because Aq is known we 
may expand 

oo 

sin0A:f(^o) = ^5j(x)P2i+incose) (A2) 

j=0 

to get the equation 

(fAAj {2j + I)i2j + 2)AAj _ H (0) R'^x sin xSj{x) 



dx"^ 



which, with the boundary condition xdAAj/dx + {2j + l)AAj = 0, has the solution 



AAj{x) 



AAj{xb)xl^+^ _ H{0)R^ 
' (4i + 3)7r2 



+ 



X 



4i+3 



„2i+l 



dx' sm x' Sj{x') ^ 1 



{x') 



M2i 



2j+l 



{X'Y^ 

dx'{x' f^+^ sin x'Sj{x') 



Differentiating Eq. (A4) and evaluating at Xh implies 
R'^ABejixb) fldAA 



X dx 



{2j + l)AAj{xb) H{d)R^xl^ r dx' sin x'Sj{x' 



+ 



where the perturbed field tangential to the boundary at Xb is 

oo 

ABe{xb,d) = Y,^Be,j{xb)P2j+i\cose) . 

3=0 



{X') 



(A3) 



(A4) 



(A5) 



(A6) 



Eq. (A5), given Eq. (A6), is used to determine the unknown parameters of the source term, 
Aj'(A). Once these parameters have been found, Eq. (A4) evaluated at x = vr determines the 
multipole moments Qj = AAj (7r)i?2-J+^: 



Qj = AA,{xby,'^' + 



4j + 3 



dx sin X Sj{x) 



Xb 



x' \ 

Xb) 



(5) 



(A7) 



Non-trivial perturbations will always engender multipole fields at the surface. 



A.l. Surface Currents from Transition Layer 

We write the zeroth order solution as Aq = Asf{x]Xb)sin^ 6, where Ag = Fq{'k) so f{-K]Xb) = 1, 
and expand 

Jmax / A \ ^ Jmax 

AJ(A) = 5^AjJ— ^^AXisinef^[fix;xb)r (A8) 
j=o ^ J=0 
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where we substituted Aq; then 



Sj{x) 



E 

n=j 



NnjAJn[fix;XbT 



27rr(n + 2)r(n + 1) 



(A9) 



r(n + j + §)r(n - J + i)r(i + i)r(-j - i) 

For sufficiently smah transition region thickness i the jump condition from Type II core to normal 
shell is given by Eq. (29), and therefore 



ABixb,e) 



TTiJb{e)H{0) sin Xb sin 61 
R 



(AlO) 



where J^b{d) is the current density parameter at the outer edge of the Type II core, which we 
compute for each field line; then with 



j=0 



(All) 



Eq. (A5) becomes 



Trejbj sinxfe i2j + l)AAj{xb) _ N^jAJn T dxs{nx[f{x; Xb)T 



R 



N. 



n=j 



{xjxbY^ 



(A12) 



which is a linear equation for the IS.Jn- The multipole moments are then found by combining Eqs. 
(ATI), (TAOI) and the solution of Eq. (|AT2|. 



A. 2. Field Lines Poking in from the Core 

We do not know A{xb^ 0) a priori. In the strong Type II core, field lines are labelled by their foot 
points xo, and our solutions determine (i) the value of 9 where they hit the boundary and (ii) the 
value of Jb{G) there, subject to whatever constraint we impose on the solution, such as for the 
"most dipolar" solution. We need A{xb.,9) in order to map jb{0) — ^coreC^); which is needed to 
determine the effect of impinging field lines on the normal shell solution. 



A. 2.1. A{xb, 



\Brl sin^e 



In this case, the background solution is Eq. (31) with current parameter Eq. (32). 
ABj{xb) = 



We take 



AAjixb). 

A key feature of this background solution is that field lines from the Type II core penetrate into 
the normal shell, bringing along the currents Jii{A) = J'b{9{A)) associated with them, where 
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d{A) = sin~\^2A/Brl). Defining AJii{A) = Ju{A) - Jshe\\{B) we have 

AJ{A) = AJu{A)e{^Brl-A) + AXhcn{A)Q{A-lBrl) 
« AJ{Ao)e{6u{x) -e) + AJ,^,ii{Ao)e{9 - Ouix)) 

where the second line assumes A ^ Aq, and 



(A13) 



Sj{x) 



[ dfiVl - fi^AJii{Ao)P2j+iHf^) 

J cos 9ii{x) 



f 

Jo 



fcos 9ii(x) 

+ / d/i^l - /x2AJshell(^0)^'2j+lH/") 



(A14) 



Generically, the first term in Eq. ( A14) is 0(1) for x — Xb^ I3eb and is 0{I3'^) for x — Xf, ~ e^. To find 
Aj7sheii(^) we must first specify its functional form and then determine the associated parameters 
from Eq. (A5) with ABgj(xb) = 0; the generic scalings suggest that Aj7sheii(^) = C'(/3). If we 
expand 

Jmax 

AJu = ^^n.ii sin^" ^ (A15) 

n=0 



on Xf), with known coefficients {AJnji} and Ai7sheii(^) as in Eq. (A8) with unknown coefficients 

{AJn,shcl\} 



Sj{x) 



1 



^ n=0 

2 



iV„,-[/3xfe/x + fix; XbT - /g"Ar„,j(x) 
(/3rfe/i? + 1)" 



AJn. 



shell 



sin^"6'ii(x) J cos eii{x) 



dMl-/i2)"+l/2p2,+/(^); 



(A16) 



although Nnj = for j > n, J\fn,j{x) / in general. The {A^„^shcii} are the solution of 



J max 



n=0 \ 

J max 

,shell 



/3"Aj7n,sheii \ r dxsmxj\fn,j{x) 



i/3rb/R + irJ J,^ {x/xbfi 

(ixsinx[/3x6/x + j{x\Xb)^ 



n=j 



{x/xbY^{Pn/R + i) 



(A17) 



A.2.2. General A{xb,i 
We can incorporate more complicated A{xb,0); let 



A{xb,e) = ^Brl 



sin^ + sin^^^ ajP2j+i"^(cos ( 



(A18) 
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the coefficients {aj} may be derived from a given A{xh,0) via inversion. Eq. (A18) implies a 
different mapping 9{A) on xi,, hence a different form of J'ii{A) for impinging field lines. With Eq. 



(A18), the background solution is multipolar. The j = component is Eq. (31) but the other 



multipoles are generically 



Aj{x) 



Brla.xl^^' 



+ 



2x2i+i ' (4j + 3)7r2 
^^■+3 n dx'sinx'sfix') 1 



^ dx'smx'sf\x') 



,2j+l 



' dx'{x' f^+^sf\x') 



(A19) 



The current density in this model is J'sheii + J^^^A), where J^^^A) must be determined. Because 
the impinging fields are weak compared with Hf,, we can determine the currents perturbatively 



via the analogues of Eqs. (A2), (A8) and (A9), plus the analogue of Eq. (A5) with the condition 
ABj{xb) = 0. 

Once the full background solution has been found, the boundary of the region to which impinging 
field lines are confined may be determined. We can then proceed as in §A.2.1 to find the perturbation 



to the background solution. Eq. (A18) implies a total flux ^(xf,, 7r/2) through the upper hemisphere. 
Consistency of the perturbation theory requires A{xi,,tt/2) « Hf^R^et- 



B. The Thin Shell Approximation for the Transition Layer 



The proton density plummets in a thin layer that extends from r = to r = ri^+l; let r = rf,+is. As 
Up falls, so does Ap and the superconductor disappears. For small enough pp^p, nonzero temperature 
will matter: for Ap ~ T, the superconductor may first become Type I before disappearing entirely 
at pf,p 7^ 0. Our treatment here assumes that T 7^ is only important in a very thin subdomain 
of the transition layer that has little effect on the magnetic field. With Ap p\ at small pp^p 
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2005 ) , Eq. ( 1 ) implies k (x pp' , so the superconductor remains 



for T = (e.g. 

Type II. However, as ki shrinks interactions among field lines become increasingly important; if 
Hci{p) = Hci{s) within the layer, then the magnetic field strength is 



(Bl) 



where Hci{s) is given by Eq. ^ and /(O) = 1 and f{z) — )• exponentially as z — )• 00. Ampere's 
law is 



ni sin XbH{0)J{A) sin dHe 



+ 



He 



dHr 



(B2) 



R ds r \ do 

if we assume that field lines do not rotate significantly, so that we can identify A with 9, and only 



retain the lowest order in i, Eq. (B2) implies 



TTis smxbH{0)J{A) 
R 



sm( 



(B3) 
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The magnetic induction field is 



B 



e dA r dA 

+ 



Ir sin 9 ds sin 9 89 ' 



(B4) 



in the thin layer approximation we assume that Br remains constant and therefore Bg = ity^l — B^/B"^. 
We assume that Hg (and Bg) remain in the —0 direction, so 



Hb 



TTis smxbH{0)J{A) 
R 



sm( 



Bl 
52 



H,i{s) + Bf 



Hcijs) 
B 



(B5) 



which is an algebraic equation for B. Once B is found from Eq. (B5), we can solve for the change 
in A: 

— ~ in sin 0^52 _ ^2 , (Bg) 



Ultimately, B approaches Eq. (B3), hence the change in A{s,9) across the layer is 

ir£smxbH{0)J{A) 



AA{9) = A{1, 9) - A(0, 9) ~ in sin^ 9 



Hb 



2R 



(B7) 



Eq. (B7) shows that AA{9) is relatively small provided that £Hb/Brb ^ 1. 
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